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Abstract 

The cosmic, general analitic solutions of the Brans-Dicke Theory for 
the flat space of homogeneous and isotropic models containing perfect, 
barotropic, fluids are seen to belong to a wider class of solutions -which 
includes cosmological models with the open and the closed spaces of the 
Friedmann-Robertson- Walker metric, as well as solutions for models 
with homogeneous but anisotropic spaces corresponding to the Bianchi- 
Type metric clasification- when all these solutions are expressed in 
terms of reduced variables. The existence of such a class lies in the 
fact that the scalar field, (j), times a function of the mean scale fac- 
tor or "volume element", a 3 = 010203, which depends on time and on 
the barotropic index of the equation of state used, can be written as a 
function of a "cosmic time" reduced in terms of another function of the 
mean scale factor depending itself again on the barotropic index but 
independent of the metrics here employed. This reduction procedure 
permites one to analyze if explicitly given anisotropic cosmological so- 
lutions "isotropize" in the course of their time evolution. For if so can 
happen, it could be claimed that there exists a subclass of solutions that 
is stable under anisotropic perturbations. 
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1 INTRODUCTION 



The first authors to realize the possibility of giving a reason for the viability of the 
"cosmological principle" without the necessity of imposing highly special initial condi- 
tions before the "inflationary programme" was developed, were Hoyle and Narlikar [|TJ. 
However, their explanation came also before the cosmic microwave background (CMB) 
radiation was discovered, and nowadays, most investigators believe that the steady state 
theory is untenable from the observations. Other investigators within Einstein's general 
relativity theory (EGR), notably Misner f2|, tried to demonstrate -unsuccesfully- that 
the large scale structure of the Universe, in particular its isotropy, could be attributed 
to the nature of the matter processes, such as dissipation, that took place at a very 
early stage of development of the Universe independent of its initial conditions (chaotic 
cosmology), that is, that the Universe lost memory of any initially imposed anisotropy 
or inhomogeneity (Barrow and Matzner [[J, Doroshkevich et. al. [|J, Misner ||, Rees 
0, Zel'dovich and Starobinsky J7|). More to the point, within theoretical cosmology 
in general and in the context of EGR in particular, one is still looking for a satisfac- 
tory explanation to the following observational facts: The formation of galaxies and 
cluster of galaxies which means that the universe is not homogeneous locally, and on 
the other hand the CMB radiation that seems to be very nearly isotropic on account 
of its Planck spectrum and the lack of structure in its intensity, from which it has been 
concluded that the large-scale structure of the actual Universe must be homogeneous 
and isotropic. The standard big bang model, thought to give the most accurate de- 
scription of the Universe, has the peculiarity that it appears to need a very special set 
of initial conditions to be viable. This state of affairs have produced several studies in 
different but related directions to obtain reasonable explanations to this conundrum, 
like the "inflationary programme" which nowadays is a popular approach, nevertheless 
not without some drawbacks, to solve also some other problems in cosmology (for a 
present review see Olive ||). A most important reason why cosmological models that 
predict inflation in the early universe are interesting is the hope they will explain the 
observed state of the universe without appeal to highly special initial conditions. Even 
so, most inflationary cosmological models have assumed Friedmann-Robertson- Walker 
(FRW) symmetry from the outset. The horizon size in the FRW models suggests the 
possibility that physical interactions could have homogenized and isotropized the uni- 
verse, and therefore that its present state could have evolved from more general initial 
conditions. Yet, it is not established if general cosmological models with non FRW 
initial conditions that actually enter, and gracefully exit, an inflationary phase initial 
homogeneity and anisotropy will be smoothed out since gravitational interactions tend 
to enhance inhomogeneities instead of smoothing them out, see Wald M. 

On experimental grounds the degree of anisotropy of our Universe is somewhat 
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bounded by the current measurements of the CMB anisotropy made by COBE ||10|| . 
Nevertheless, present observations are still far from elucidating all the properties that an 
actual model of the universe should have. Even the degree of anisotropy of the primeval 
radiation or the counts of radio sourcess or galaxies in the various directions in the sky 
involve some uncertainties and, agreement between the observed chemical composition 
and the predictions of the Friedmann models merely signifies that the time rate of 
change of the volume occupied by matter and the rate of expansion of a Friedmann 
model concide. On the other hand, astronomical observations seem to imply that 
isotropization must occur at a rather early epoch, maybe even for z » 10 9 . It is feasible 
that for particular anisotropic models this could be achieved for the instantaneous values 
of the anisotropy parameters. But, the case is that there are observational properties 
that depend on the degree of anisotropy over an extended period of time, again like the 
CMB isotropy. 

In EGR, without the aid of a cosmological constant or inflation, Collins and Hawking 
|jlT[1 examined the question in terms of an "initial conditions" analysis. They obtained 
that the set of spatially homogeneous cosmological models approaching isotropy in the 
limit of infinite times is of zero measure in the space of all spatially homogeneous models, 
which in turn implies that the isotropy of the models is unstable to homogeneous and 
anisotropic perturbations. Yet, other authors define the concept of isotropization in 
a different way (see Novikov [JT^] , Mac Callum |T3| and Zel'dovich and Novikov \T%\). 
Therefore, in the literature concerned with the mathematical analysis of anisotropic 
models the term isotropization is often mentioned but its precise definition is author 
dependent. For the Bianchi-Type models it is claimed that a positive cosmological 
constant provides an effective mean of isotropizing homogeneous universes (Wald ||). 
However, in this context Barrow has shown that contrary to previous expectations, 
perfect fluid cosmologies need not approach isotropy and homogeneity as t — ► oo. 

Following the general line of thought on this subject put forward by Zel'dovich and 
Novikov, one declares that homogeneous cosmological models that isotropize are those 
that "approach to a Friedmann model in the course of time as the universe expands" , 
which means that "its geometric and dynamic parameters as well as those concern- 
ing the distribution and motion of matter and radiation are nearly the corresponding 
quantities in a Friedmann model" . Accordingly, one assumes that the idea conveyed by 
isotropization is the property that an arbitrary solution possesses to model our Universe 
which permits it to evolve from an initial, general state, into a state that is presently 
isotropic on a "large scale" and is, therefore, well described by a Friedmann solution. 
We reduce the general scope of the problem by assuming initial homogeneity, limit- 
ing ourselves to test the isotropization properties of certain specific non inflationary 
solutions in the Jordan-Brans-Dicke (JBD) cosmological theory. So, if at its outset 
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the universe was not in isotropic expansion, the above ideas imply that one can exam- 
ine, on a first approximation, the properties of homogeneous but anisotropic models 
assumed to describe correctly the early stages of its expansion. Of these, only those 
Bianchi-Type models whose group type comprise FRW models may isotropize: Types 
I, V, VII , Vll h and IX. 

In this paper the above concept of isotropization is dealt with in a direct, but 
admittedly limited, way by qualifying and quantifying it through a " Raychaudhuri 
type" equation common to all Bianchi-Type models: Given an explicit solution, one 
can directly check if it may or may not approach to a Friedmann regime in the course of 
its cosmological time evolution, specifically, if the different anisotropic scale factors of 
a Bianchi model in the various directions approach arbitrarily near to a unique, single 
function of time. By this procedure one can then answer the question, at least for 
some representative spatially homogeneous models of the Bianchi- Type (I, V, IX), of 
whether, and if so, how in the JBD cosmological theory a present large scale isotropy 
resulted from an initially anisotropic but homogeneous expanding universe. 

2 FRW FIELD EQUATIONS 

The JBD field equations for the FRW cosmology with a barotropic, perfect fluid, 
p = (3p , — 1 < (3 < 1 (the (3 — |, equation of state for incoherent radiation or 
ultrarelativistic matter is excluded) are 

po 3(1+/3) = M/3, Mp = const 

The dynamic equation is 

[2(2 - 3(3) + 3(1 - (3) 2 oj\mp -6(1 - (3)k 

^ ~ a^-W ' 

and the constraint equation is 

3 (r\ _ i /v\ 2 _ (i - 3/3) / (l- 3/3)7^77+^ /y\ 
2(1-/3) (i-/?) 2 Wy (i -/5) 2 v i> )W) 

[2-3/3 + f(l-/3) 2 io] / (l- 3/3)777.^ + 77 \ 2 3[2 + io(l-/3)(l + 3[3)}m p ( 4 ) 
(1-/3) 2 ^ i, ) + 2(1-/3)^ 

= , 



8ttM, 



/3 



3 + 2uj 



(1) 
(2) 



k = 0, ±1 (3) 
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where i/j = 0a 3 ( 1- ^, is the JBD scalar field, a is the scale factor, ui the coupling 
parameter of the theory, 77 the "cosmic time parameter" , 770 an integration constant, and 
( )' = d v , where dt = a 3/3 di] (for details see Chauvet and Pimentel [Hi], and references 
therein) . 

For k = equation (|3|) is directly integrated. One gets 

ip = Arj 2 + Br] + C (5) 

where A, B and C are constants such that 

3 



.4 



m p . (6) 



Substitution of (§) and (||) in the constraint equation hands out the following results. 
The constant B is undetermined and so, up to B, C also remains undetermined. There- 
fore, three different possible cosmic solutions to the FRW flat space (k = 0) exist dis- 
tinguished by the sign of the determinant, A = B 2 — AAC, which itself depends on 
the relation between the equation of state, through ft, and the coupling parameter u in 
a rather complex way. The behavior of the scalar field implies, for each type of deter- 
minant, A > 0, A < and A = 0, the possible existence of two branches: essentially, 
ones with an increasing function of time and the others with a decreasing function 
of time (the solutions are given explicitly and thoroughly discussed by Gurevich et.al. 



17], Ruban and Finkelstein |18| and Morganstern [19 



is obtained by the straightforward integration of 

(f/_ _ (1 - 3/3)77^77 + T] 

1p 

and the scale factor, a , found from it through the definition of ip: 



(7) 



a 3(l-/3) 



Arj 2 + Bri + C 



Equation (|j) does not involve the curvature constant k explicitly and so, ip = 
Arj 2 + Br/ + C is also a solution to actually both the open and closed space dynamic 
equation ([3D provided that 

= 2 + (l-.3)(l + 3/?V 
y 2(l + 3/?)Jfe p V ; 

The same as in the flat space case, A, B and C are obtained from the constriction 
equation ((4[). For both, k = +1 and —1, it is valid that 

-(1 - 3/3) 2 m /3 
(1 + 3(3) ' 
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B = - 2 (lri)*' ( 10 ) 

c = 711 



Its determinant is then 

A = B 2 -AAC = 0. (11) 
The explicit solutions for these two models are 

-1 f[2+ 1-/3 1 + 3MV^ r/i om , i-^ 
' j [(l-3/3)m^+77o] i-35, 

(12) 



r (1 + 3/3)777^ \ -2jfe 
and 



-2k \ 2 ( 1 " 3 ' 3 ) i 

[(1-3/3)7^ + 770]^ . (13) 



V[2 + (l-/3)(l + 3/3V]m2; 
The solutions for the JBD flat space were previously obtained by several authors 



(Gurevich et. al. \Tf\, Morganstern [19) and references therein, and in another context 
by Chauvet and Pimentel [16|). 

Next, we present the anisotropic Bianchi field equations in the above variables in 
order to analyze later their asymptotic solutions. 

3 ANISOTROPIC FIELD EQUATIONS 

Three extra equations, and simple modifications to the FRW equations (0)-(§) pre- 
sented above, describe the Bianchi-Types I, V and IX examined in this paper. Equa- 
tions and (0) remain formally the same, while equation (||) gets its "curvature" term 
modified and is then written 

w\ _ [2(2-3^ + 3(1-/^ = (1 _ ^ , R> (14) 

The constriction equation is a " Raychaudhuri type" equation so that the left hand side 
of equation (|) remains unaltered, but instead of being equal to zero as in the FRW 
cosmology, it is in this case: 

3 /y\ 1 fij'V (1-3/3) f(l-3p)mpr) + r) \ (V\ 



2(1-/3) J (1-/3) 2 V^/ (1 — /3) 2 V V> ) \i> ) 

[2-3/3 + |u;(l-/3) 2 ] / (l - 3/3)7^7/ + Vo \ 2 3[2 + u(l - /3)(1 + 3/3)]m p (15) 
(I-/?) 2 V i> ) 2(1-/3)7/; 

= -(#! - H 2 f - (H 2 - H 3 ) 2 - (H 3 - H,) 2 = a{rj) . 
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a is, for short, the "shear", a = 0, is a necessary condition to obtain a FRW 



cosmology since it implies Hi = H 2 = H% ( see Chauvet et.al. [20|). If the sum of 
the squared differences of the Hubble expansion rates tends to zero it would mean that 
anisotropic scale factors tend to a single function of time which is, presumably, the 
scale factor of a corresponding Friedmann model. However, in general not all Bianchi 
models contain a FRW space-time. 

The three extra equations describe the dynamical evolution of the "anisotropic scale 
factors" ai,a,2 and a 3 : 

(if) Hi)' = [1 + {1 - P)u]mp + ip a 60 *Rij . i = 1, 2, 3 (16) 

From equations ([14]) through (fL6"|), and for the rest of this paper, we use the following 
notation and conventions: a is presently the mean scale factor, a 3 = a^a^, the H^s 
i = 1,2,3 are the Hubble expansion rates, Hi = d i ja^ *Rj is the "spatial three- 

3 

curvature" that belongs to a given Bianchi-Type model, *Rj = is a column 

sum, and the *Rij are "partial curvature" terms pertaining to specific scale factor 
dynamic equations, in our case: 

IV IX 

2/a\ [a\ -a A 2 -a\ + 2a 2 2 aj] / (-2a 6 ) 

(17) 

*R ij = 2/a\ \a\ - a\ - a\ + 2aja 2 3 }/ (-2a 6 ) 

2ja\ \a\ -a\-a A 2 + 2a\a\\l(-2a 6 ) 

In the next section it is first shown that ip = Arf + Bt] + C includes solutions for the 
above, purposely chosen homogeneous but anisotropic, cosmological models. It will be 
shown that the obtained ip class of solutions consist of two parts: the isotropic and the 
anisotropic one. Then, it will be clear that the latter approach to zero as the cosmic 
time parameter evolves, i.e., these solutions tend asymptotically to their corresponding 
isotropic group solutions, the FRW models. 

4 ANISOTROPIC SOLUTIONS AND THEIR 
ASYMPTOTIC BEHAVIOR 

We show next that ip = Ar] 2 +Br]+C is a solution for the, homogeneous but anisotropic, 
Bianchi-Type cosmological models. 

Bianchi- type I 
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For this Bianchi type model, *Ri = 0. It is direct to see by substituting 



ij) = Arf + Brj + C 



into equation (H), that A has the same expresion as the one given by equation 
(from now on we attach a subindex to the A, B and C to distinguish between the 
different Bianchi models): 

Ai = [2 -3/3+ |(1 -P) 2 u}mp ■ (18) 

By direct substitution of the above results into equation ( |15|) one finds that B l 
remains undetermined and may be put equal to any convenient, but arbitrary value, 
and that 

3(1 - (3) 2 (3 + 2u)m p C I = 



(hi + hZ + hi)-m TiZA I + B/ + (l-3p)m l3 Ti o B 1 , (19) 



where the h^s are constants such that 

1 a 3 ' hi , , 

The non-vanishing constants hi's determine the anisotropic character of the solutions. 
They obey the condition 

hi + h 2 + h 3 = 0. (21) 

By integration of equation (|20|), using equation (0) and equation (|3p, one finds 
explicitly a = a(rj) (first obtained by Ruban and Finkelstein [IS|, see also Chauvet and 
Guzman [BTJ and Chauvet 



For the Bianchi's in general, when equation (||) is substituted into the Raychaudhuri 
equation (15), one obtains the shear as a function of ip and the h^s, 



KU 2^ 2 v ; 

This term permits, beside the only one allowed A = solution for the k ^ 0, FRW 
models, two other solutions with A^O such that a — > as i] — > oo (or t — ► oo). It is in 
this sense that these solutions may isotropize in the course of their time evolution (note 
a mathematical characteristic of the anisotropic solutions shown by the above results, 
i.e., the relation between the exponents hi> s and the B and C coeficients). 

Bianchi- type V 
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Equation (|5|) is a solution for this Bianchi model, with A, B and C equal to 

(1 - 3/3) 2 m /3 



A 



-2 



;i + 3/5) 

1 + 3/3 J 



(23) 



*7o 



A v and I? v are equal to the ones obtained for the isotropic, k = ±1 cases, but 



So that, 



m p (1 + 3/?) C v 



[2 + (l-/3)(l + 3/?)a;]m /3 



;i + 3/3) 2 (fri + frj + /^) _ 2 
18/3 + u;(l + 3/?) 2 7/0 



(24) 



and 



-2(1 + 3/5) 



1 a 3 ' a'i 



3(1-/3) 
2(1-3/3) 



A^ 2 + S F?7 + C v 



^1 - 3f3)mpr] + T] 



1+3/3 
"2(1-3/3) 



3a3 a x (l + 3(3)(A v ^ + B vV + C v ) 



The scale factors are 



-2(1 + 3/5) 



2(1-3,3) 



A^ 2 + B vV + C v 



and 



or 



a 2 = a,i exp 



a 2 = ai exp 



[2 + (l-/3)(l + 3/5V]m^ 

-2/i, . f-2(l - 3y9)[(l - 3P)m p rj + r]o} 



2(1-3/3) 



(25) 



(26) 



(27) 



A 

2fr 2 



•arctanh 



arctan 



(l + 3/3)VA 

-2(l-3/3)[(l -3^77 + 7/0] 
(l + 3/3)>/=A 



with 



CL2CL3 — Q-i 

These solutions are new. However, for A = 0, 

CI3 = Q2 ~ ^1 



, A > (28) 

,A<0 (29) 
(30) 

(31) 



is obtained. This last solution is clearly seen to be the one previously obtained for the 
isotropic, FRW model, with an open space {k = —1). Again, the A = is obtained 
only if hi — hi — h% — 0. 

Independent of the value A might have, h\ + hi + /13 = is always true. In the 
Type V models with A ^ one must have that h 2 = —h 3 with hi = 0. For the latter 
case, truly anisotropic solutions are obtained with 

-8(1 - 3/5) 2 



A = B y - AA V C V 



18/3 + (1 + 3p) 2 u 



2, ,2 



hi . 



(32) 
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C v , being proportional to the sum of the squares of the constants h 2 and h 3 , carries 
the information concerning the nature of the anisotropic character of this Bianchi-Type 
model. 

Since equation (|22| ) holds for all the Type V models, the A ^ solutions could 
have had asymptotic behaviors to call them "nearly isotropic in appearance" if cr — > 0, 
when 7] — > oo. In this regard one finds that: for the A < solution a 2 tends to 
axe.x'pijshil V '—A) and 03 tends to a\exp(— Kh<ij\J — A) so, when r\ reaches the value 
(100v^— A — B V )/2A V these two scale factors differ from each other by one percent, 
and this solution is then "ninety nine percent near" the A = solution, corresponding 
to the FRW cosmology for k — — 1 (see equations (0) and (0))- On the other hand, 
for the A > model the scale factors can never approach to a same, single, function 
of rj. The reason for this is that r\ is bounded, see equation (p8|). Nevertheless, it is 
significant that this last can be an inflationary solution. 

Bianchi type IX 

ijj given by equation (g) is, likewise, a solution for this model with the H^s given 
by equation (p0|). However, in this case the h^s cannot be constants. Instead, the 
h[s = hi(rj)'s are now new, and unknown, functions of r\. 

With equation (|J) substituted into 

[Hii/i]' = [1 + (1 - p)u]mf, + a 6/ V *R iIX , i = 1, 2, 3 (33) 

one must solve for 

, R 2A TX - [2(2 - 3/3) + 3(1 - (3) 2 u]m n , 
K = a e ^*R iIX + IX 3(1-/3) ■ 1 = 1,2,3 (34) 

The sum of the above three equations, 

a«V*,„ = U " - P(2 -f +; 3(1 ' (35) 

(1 P) 

is, given explicitly in terms of ai, 02, and (13, 

a\ + a\ + a\- 2{a\a\ + a\a\ + a\a$ _ [2(2 - 3/3) + 3(1 - (5) 2 uj]m p - 2A 



IX 



2a 6 ( 1 -« (1 - 

(36) 

from which any chosen scale factor can be solved as function of the other two remaining 
ones. 

On the other hand, equation ([15]) gives 

h\ + h 2 2 + h\ = K 2 = ~ 2(1 ^ )2 [^V + 0*7 + g] , (37) 
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where the P, Q and S constants, given in terms of A IX , B IX and C IX , stand for 

P = XA IX -[AA IX -Y](l-3(3) 2 m p , (38 



and 



where 



and 



Q = XB IX - [AA ix t]q - 2Ym p r ]G + 2(1 - 3/3)5 J (1 - 3(3)m p , (39) 

S = XC IX - [2A + 2(1 - 3(3) mf3Vo B IX - Ym 2 V 2 } ; (40) 

X = 3(1 + 3/3) (1- /jyWi/j + 6(1 - P)m p - 2(1 + 3(3)A IX , (41) 

y = 2(2 - 3/3) + 3(1 - f3) 2 u . (42) 



The isotropic model solution that belongs to this Bianchi-Type model is obtained 
when A = 0, where one has that hi(j]) = h 2 (ri) = /i 3 (77) = (see section 2). 

Under any circumstance the functions h^s, which must still obey the condition 
hi + ti2 + h% = 0, determine the anisotropic character of the solutions. 

The h^s can be given as 



K 2 + AK + 1 



and 



3{k 2 + k + 1; 

-k 2 + 2k + 2~ 
3(k 2 + k + 1)_ 

2k 2 + 2k - 1" 



K 



K 



K 



(43) 



(44) 



(45) 



_3(k 2 + k + 1)_ 

where now k is another, new and yet unknown, function of 77: Unfortunatly, for A ^ 
we were not able to obtain the explicit functional dependence of k = Even so, an 

asymptotic isotropic behavior, similar to other models, for the present solutions is also 
expected based on the strength of equations ( [22]) and (|37|). 



5 DISCUSSION AND CONCLUSIONS 

The JBD cosmological equations for perfect fluids with barotropic equations of state 
is seen capable of beeing displayed, through the use of reduced variables, in a way 
which first permits one to obtain non trivial, significant solutions with little effort 
and next, but more important, to express them in terms of the single function if) = 
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Arf + Brj + C. The fact is that the aforementioned solutions belong to a class which 
embraces Bianchi-Type models some of which, in turn, comprise the FRW isotropy 
groups. Moreover, stated explicitly this class contains the general (analytic) matter 
solutions for the Bianchi-Type I model as well as solutions for the other two Bianchi- 
Types examined in this paper, which in turn include special ones, a subclass, that 
tend aymptotically, as rj — > oo, to corresponding FRW solutions. The reason for the 
existence of this set is that the functional form of the product of the scalar field, 0, times 
a power of the mean scale factor, a 3 = a^a^, as a function of the time parameter r] is a 
solution to the equations used in this work independent of the metrics that give rise to 
any possible present anisotropy for Bianchi-Types I, V and IX models, and it has the 
FRW form. In other words, for a perfect fluid with a barotropic equation of state, we 
have shown that there exists a class of solutions for the Bianchi I, V and IX types that 
contain their corresponding FRW models. The Type V solutions are new, as well as 
those for the Type IX, but in the latter case because of the complexity of the curvature 
terms it is only possible to give the explicit form, in terms of r], of the scale factors a^s 
up to the single, unknown, function k = Nevertheless if 77 — > 00, an asymptotic 

isotropic behavior for the A ^0 solution should be expected in view of equations (p2"D 
and (|3?D . Moreover, there are also other solutions obtained from ip, through equation 



), which describe other Bianchi-Type models, see Chauvet and Guzman [21]. All of 



the above are salient and remarkable properties of the class of solutions that we have 
found. Even more, the ip solutions also include the asymptotic solutions, these are the 
Type I models (A = 0, > and < ), for all Bianchi-Types near the initial singularity, 
when the spatial curvature terms can be neglected. They are in this sense -up to the 
matter terms- comparable to the Kasner vacuum solution of EGR. 

At this point we want to mention that, in JBD, the separate cosmological models 
are to be distinguished between themselves in a first instance, through the different 
values that the constants A, B and C obtain. Also significant is the fact that B and 
C carry the physical information on the nature of the presence, or even the absence, 
of the anisotropy that any given models may have. Remind that these constants also 
determine the value of the discriminant A in terms of the physical parameters (3 and 
to. We stress the fact that solutions for the Bianchi models corresponding to A = 
that possess an isotropy group, recover the FRW solutions: Type I goes into the 
corresponding flat FRW one, Type V into the open FRW one and the Type IX into 
the closed FRW one. Meanwhile, for the anisotropic solutions with A / the models 
isotropize in the course of their time evolution by tending in the Type I case to its 
corresponding FRW model, while for Types V and IX their evolution is toward the 
A = solution which is the only one available to describe a non-flat FRW cosmology 
in the JBD context. 
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The present observational evidence points to a high degree of isotropy of the, as- 
sumed relic, CMB radiation and, if so, it is a decisive argument in favour of the, nowa- 
days, large scale homogeneous and isotropic expansion of the Universe. Then, if the 
initial stages of the expansion had a homogeneous but anisotropic behavior one could 
follow, within the JDB, how an actual nearly isotropic expansion can come about. 

Acknowledgment This work was supported by CONACYT grant 400200-5-3672E. 
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